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ABSTRACT 

We present a new study of the collisional runaway scenario to form an intermediate- 
mass black hole (IMBH, Mbh > 100 M 0 ) at the centre of a young, compact stellar 
cluster. The first phase is the formation of a very dense central core of massive stars 
(M* ~ 30 — 120 M 0 ) through mass segregation and gravothermal collapse. Previous 
work established the conditions for this to happen before the massive stars evolve off 
the main sequence (MS). In this and a companion paper, we investigate the next stage 
by implementing direct collisions between stars. Using a Monte Carlo stellar dynamics 
code, we follow the core collapse and subsequent collisional phase in more than 100 
models with varying cluster mass, size, and initial concentration. Collisions are treated 
either as ideal, “sticky-sphere” mergers or using realistic prescriptions derived from 
3-D hydrodynamics computations. In all cases for which the core collapse happens 
in less than the MS lifetime of massive stars (~ 3Myr), we obtain the growth of a 
single very massive star (VMS, M* ~ 400 — 4000 M 0 ) through a runaway sequence 
of mergers. Mass loss from collisions, even for velocity dispersions as high as a v ~ 
lOOOkms -1 , does not prevent the runaway. The region of cluster parameter space 
leading to runaway is even more extended than predicted in previous work because, 
in clusters with a v > 300kms _ , collisions accelerate (and, in extreme cases, drive) 
core collapse. Although the VMS grows rapidly to > 1000 M 0 in models exhibiting 
runaway, we cannot predict accurately its final mass. This is because the termination 
of the runaway process must eventually be determined by a complex interplay between 
stellar dynamics, hydrodynamics, and the stellar evolution of the VMS. In the vast 
majority of cases, we find that the time between successive collisions becomes much 
shorter than the thermal timescale of the VMS. Therefore, our assumption that all 
stars return quickly to the MS after a collision must eventually break down for the 
runaway product, and the stellar evolution of the VMS becomes very uncertain. For 
the same reason, the final fate of the VMS, including its possible collapse to an IMBH, 
remains unclear. 

Key words: Galaxies: Nuclei — Galaxies: Starburst — Galaxies: Star Clusters — 
Methods: N-Body Simulations, Stellar Dynamics — Stars: Formation 


1 INTRODUCTION 

In our first paper (Freitag. Rasio. fc Baumgardtll200fil : here¬ 
after Paper I), we have presented our numerical approach 
to the study of stellar collisions in young, dense star clus¬ 
ters with a broad stellar mass spectrum. It is based on 
the use of ME(SSY)**2, a Monte Carlo code to simu- 
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late the long-term evolution of spherical clusters subject 
to relaxation, collisions, stellar evolution and taking into 
account the possible presen ce of a central massive object 
(Freitag fc Benzll200ll . f2002i) . Our main motivation is to in¬ 
vestigate in detail the runaway growth of a massive object 
(“very massive star”, h ereafter VMS) during core collapse 
lEbisu££dc^^iQj20(IlF[lYrte! ri es Zwart fc McMillanl Eo02 i: 

lffl^iuF^Veita^fc^siol2n(M 'hereafter GFR04). This pro¬ 

vides a natural path to the formation of an intermediate- 
mass black hole (IMBH) in a dense star cluster, or a seed 
black hole (BH) in a proto-galactic nucleus In Paper I, 
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we presented a number of test calculations to validate our 
Monte Carlo code and compare its results for simple ide¬ 
alized systems to those from A-body, Fokker-Planck and 
gaseous-model codes. These comparisons established that we 
can reliably follow the relaxation-driven evolution of clusters 
with a broad mass function, all the way to core collapse. We 
also found good agreement with models of dense clusters in 
which mergers between stars create a mass spectrum start¬ 
ing from a single-mass population, th us accelerating col ¬ 
l apse and leading to collisional runaway llQuinlan fc Shapirol 


Here we use ME(SSY)**2 to perform a large number of 
calculations covering broadly the parameter space of young 
star clusters. Specifically, we vary systematically the num¬ 
ber of stars (in the range 10 5 — 10 s , using up to 9 x 10 6 
particles) and the size of the cluster (with half-mass radius 
values from 0.02 to 5pc), and we consider systems with low 
and high initial concentration (King parameter Wo = 3 and 
Wo = 8). Our work is guided by the main finding of GFR04, 
namely that for clusters with realistic initial mass functions 
(for which the ratio of the maximum stellar mass to the av¬ 
erage mass is larger than ~ 100), mass segregation drives the 
cluster to core collapse in a time not longer than 15% of the 
initial central relaxation time. It is therefore expected that 
in any cluster for which this relaxational core-collapse time 
is shorter that the main-sequence (MS) lifetime of the most 
massive stars (M* fts 100Mg), i.e., ~ 3Myr, these objects 
will eventually collide with each other in the liigh-density 
core, likely in a runaway fashion. As we show in Sec. |3 our 
results confirm this expectation. 

Our work on the subject was inspired by the in¬ 
vestigations of Portegies Zwart and collaborators who 
revived the study of the runaway scenario t hanks 
to A—body s imulations ijPgrteeies Zwart et a.lJ ll 9991: 
IPorteeles Zwart Rr. McMilla.nll2002Fl . Although many key as¬ 
pects of the process were already known from older works 
(in p articular ISpitzer fe Sasla\Jl966l:IColga,tdll967tlSanderJ 
p97a iBeeelman fc Rees OH; IVishniad ITH [H Eli! 
[Quinlan fc Shapirol Il99( ; see paper I for a review of the 
field), these recent A—body studies were seminal in con¬ 
sidering the effects of collisions in a cluster with a real¬ 
istically broad IMF. This demonstrated explicitly for the 
first time how mass segregation can lead to a collisional 
phase by concentrating massive stars in a small central vol¬ 
ume and made it clear that the post-MS evolution of the 
most massive stars initially present can prevent the runaway 
phase by driving cluster expansion. Furthermore, Portegies 
Zwart and collaborators showed that, contrary to the some¬ 
what unrealistic situation in single-mass clusters llLedll987t 
iQuinlan fc Shapirolll99oll . dynamically formed binaries, far 
from preventing the runaway phase by halting the central 
concentration increase, promote it, as many binary interac¬ 
tions lead to stellar mergers. As noted in paper I and sug- 
gested by A—body simulatio ns with higher particle numbers 
dPortegies Zwart et ali t2flQ4 ), in a broad-IMF cluster with 
IV* > 10 6 , binaries probably cannot form through 3-body 
interactions before the collisional phase is reached. 

It is important to stress, however, that our MC simu¬ 
lations probe a different regime from that explored by the 
direct A—body approach. We consider systems with higher 
number of stars in the central region (inside a few core 
radii). For this reason, some of our most important results 


are only superficially similar to the findings of Portegies 
Zwart et al. In particular, in such large-A clusters with a 
long initial central collision time, a genuine core collapse, 
itself driven by mass segregation, appears to be the condi¬ 
tion for collisional runaway. The situation for systems with 
a smaller number of stars is more complex because there 
may not be enough massive stars in the central regions 
to d rive suc h a core coll apse jPorte gies Zwart et alJ 
^2^_Portegies Zwart fc McMillan 2Q22,; 


MbMilian fc Portegie^Zwa^TT^OM : IPo^gles Zwart et al'l 


2004 . 


Our paper is organised as follows. In Section 2, we 
present the simulations we have performed and explain their 
results. In Section 3, we summarise our findings and discuss 
avenues for future research in the field. Detailed presenta¬ 
tions of the runaway scenario, the physics at play, our nu¬ 
merical methods and test computations are to be found in 
Paper I. 


2 SIMULATIONS 

2.1 Initial conditions and units 

In this work, we consider the evolution of isolated spherical 
stellar cluster with a broad, realistic initial mass function 
(IMF), subject to two-body relaxation and stellar collisions. 
When not stated otherwise, we use a Salpeter IMF, for which 
the number, dA*, of stars with masses between M* and M , + 
dM * is given by 

^ oc M~ a for M min < M„ < A/ max (1) 

dM* 

with a = 2.35, M m i n = 0.2 Mq and M max = 120Mg. For 
this stellar population, the average mass is (M*) ~ 0.69 Mq. 
There is no initial mass segregation. To investigate the role 
of the initial concentration of th e cluster, we consider King 
models with Wo = 3 or Wo = 8 dBinnev fc Tremainelll987I 
Sec. 4.4). We do not enforce tidal truncation because it was 
shown in GFR04 that this does not affect the central regions 
which completely dominate the evolution. 

We refer to GFR04 (Section 3 and Table 1) for detailed 
explanations about the important physical parameters of 
such clusters and units. In keep with the tradition, when 
not stated ot herwise, we are using the A-body unit system 
(lHenonl !! 197lf) defined by G = 1, M c i(0) = 1 (initial total 
cluster mass) and C/ c i(0) = —1/2 (initial cluster potential 
energy). As time unit, we prefer the “Fokker-Planck” time 
Tfp to the A-body unit Tnb because the former is a re¬ 
laxation time while the latter is a dynamical time; they are 
related to each other by Tfp = A*/ln(7 c A*) Tnb where A* 
is the total number of stars in the cluster. As explained in 
paper I, we conservatively use 7 C = 0.01 to determine re¬ 
laxation times in this work. For King models, the A-body 
length unit is close to the half-mass radius, Anb — 1.19 Rh 
for Wo = 3 and Anb — 1.15 Ah for Wo = 8. The core radius 
(Spitzer87, Eq. 1-34 and Paper I) is A CO re — 0.543 Anb and 
encompasses a fraction 0.238 of the total mass for Wo = 3; 
for Wo = 8, these values are 0.121 Anb and 0.0531, respec¬ 
tively. 

Expressions for half-mass and local relaxation times are 
given in GFR04 and Paper I. We denote the initial values 
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Figure 2. Time of collisional runaway occurrence as a function 
of central initial collision time for 120 Mq stars (Eq. |3), for all 
simulations with standard Salpeter MF where runaway happened. 
Times are given in units of the initial central relaxation time. 
Triangles and circles are for Wo = 3, 8 models respectively. The 
squares are Wo = 9, 12 models for the cluster MGG-9 in M 82 (see 
Sec. 12.5.^1 . These simulations were carried out with a Kroupa MF 
extending from 0.1 to 100 Mq. Solid symbols correspond to the 
“sticky-sphere” treatment of collisions, open symbols to the SPH 
prescriptions for collision outcomes. The dashed line is an ad-hoc 
parameterisation of the sticky-sphere results given in Eq.|l] Note 
that for large collision times, one gets t ra — 0.12t rc (0), in good 
agreement with GFR04 for the relaxational core-collapse time 
but with a large dispersion in the result. The star corresponds 
to a simulation with no 2-body relaxation but the same initial 
conditions as those for the triangular point below it (see Fig. EJ. 


of the half-mass and central relaxation times by f r h(0) and 
f r c(0), respectively. 

Table 1 lists the initial conditions for all runs performed 
in this study. N p denotes the number of particles used in the 
simulation. We also give in this table the value t ra of the time 
when runaway started for all runs in which a VMS formed. 


2.2 Overview of the results of the standard set of 
simulations 

We first establish the conditions for runaway formation of 
a VMS using the simplest and most favourable prescription 
for the outcome of collisions, i.e., assuming that every time 
two stars come closer to each other than the sum of their 
radii, they merge without any mass loss. With such a clean 
collision physics, one expects the runaway to happen pro¬ 
vided that: 

(i) Core collapse, driven by mass segregation due to 2- 
body relaxation (as studied in GFR04), occurs within t * ~ 
3 Myr or the cluster is initially sufficiently collisional. 

(ii) The collisional cross section is a steeply increasing 
function of a star’s mass, S co ii oc M* with r/ > 1. 


What “sufficiently collisional” means is not easy to estab¬ 
lish from simple principles. One may think that a sufficient 
condition is that the collision time at the centre, for a star 
near the top of the IMF (i.e., 120 Mq for our standard IMF), 
tcoiiljvf (see Eq. [3J, be shorter than the MS life-time for 
such a star. But this does not account for the structural evo¬ 
lution of the cluster due to relaxation and collisions them¬ 
selves. A nother condition is suggested by the coagula tion 
equation (iLedll993l . I200CX -: . Malvshkin fc Goodmanll200ll : see 
also the simplified mathematical analysis encompassed in 
equations 6 and 7 of Paper I). In the case of strong gravita¬ 
tional focusing, it reduces to R* oc M% with /3 > 0 ( R * is the 
stellar radius), a condition always obeyed by MS stars (but 
not during the pre-MS stage). If gravitational focusing is 
strong (i.e., for velocity dispersion a v > 500— lOOOkms -1 ), 
one needs /3 > 0.5, which is (although marginally) not sat¬ 
isfied by our standard VMS M-R relation (/3 = 0.47, see 
Paper I). In any case, the relevance of the second condition 
is not clear because it stems from a type of analysis which 
disregards the stellar dynamics, most noticeably the role of 
mass segregation. 

Fig. □ represents the (M c i, .Rnb) plane in which are 
plotted, for Wo = 3 and Wo = 8, the conditions t C c| rlx = 
0.12f rc (0) = t* = 3Myr and tcoii| Mmax = 3Myr. fcc| rlx 
denotes the relaxational core-collapse time, when collisions 
are absent. The value t C c| rlx = 0.12t rc (0) is the average 
value yielded by ME(SSY)**2 when collisions are absent 
(see Paper I and discussion below). We also indicate the 
initial one-dimensional velocity dispersion at the centre, 


<T 0 = cr„ jC (0) ~ 0.55<J 


~ 36 kms 


GM C i 
Rnb 

- i ( M c i 


V 10 6 M 0 


1/2 



( 2 ) 


This is approximately correct for any Wo value. By com¬ 
paring (To with the escape velocity from a stellar surface, 
500 — lOOOkms - , one can estimate whether gravitational 
focusing (included in all simulations) and collisional mass 
loss (not accounted for in the sticky-sphere approximation) 
may play an important role. 

From this diagram, it is clear that only relatively small 
clusters may have t cc | rlx < 3 Myr without being initially 
collisional (t co ii| M > 3Myr); for this one must have 
N * < 10 6 for Wo = 3 and N* < 3 x 10 6 for Wo = 8. This fact 
was not included in GFR04, who considered only the fcc| rix 
condition. Consequently, there are in principle more clusters 
which may form VMSs through collisions. Even though ini¬ 
tial conditions with collision time shorter than a few million 
years are questionable, they are an interesting idealisation to 
pave the way to more realistic simulations of the collisional 
formation of dense clusters. 

We show in Fig. 0 the initial conditions for all sim¬ 
ulations carried out with the standard IMF and collisions 
treated as pure mergers or with the SPH-generated prescrip¬ 
tions. We also indicate the outcome of each simulation, i.e., 
whether a VMS grew or whether the core collapse was ter¬ 
minated by stellar evolution of the massive MS stars before 
this could happen. We set the lower limit for successful VMS 
growth at My ms Js 400 Mq . This is a relatively arbitrary 
value which matters only in the case of the four simulations 
with Wo = 3, Nt, = 10 5 and Rnb = 0.5 pc. These runs were 
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Number of stars 



Total cluster mass [M 0 ] 


Figure 1. Diagram summarising the initial conditions and outcomes of the “standard” cluster simulations with collisions treated as 
perfect mergers or using the SPH prescription. Each point represents the cluster mass and (initial) AT-body length unit (%b) for one 
simulation. Triangles correspond to Wo = 3 King models, round dots to Wo = 8. The standard 0.2 — 120 Mq Salpeter IMF was used 
in all cases. The solid diagonal lines (of negative slope) show the condition for the relaxation-driven core collapse time to be 3Myr, 
£ cc | rlx = 0.12t rc (0) = 3Myr for Wo = 8 (top) and Wo = 3 (bottom). The long-dashed lines (of positive slope) show where the central 
collision time for a 120 Mq star (with any other star) is 3 Myr. The dot-dashed curves indicate were the time for occurrence of runaway 
should be 3 Myr, according to our parametrisation (see Ea. IH and text). Short-dashed lines indicate approximately the 1-D central 
velocity dispersion. Solid symbols are for simulations which resulted in runaway formation of a VMS; open symbols are for cases in which 
stellar evolution interrupted core collapse before a VMS could grow. In some cases, several simulations with the same 1 ?nb and A* (but 
different Wo, initial realisations of the cluster, random sequences for the MC algorithm, number of particles or collision prescription) 
were carried out. These cases correspond to symbols connected together with a thin line. Asterisks indicate simulations done with a 
number of particles smaller than the number of stars (3 x 10 5 or 10 6 particles for Wo = 3 cases, 3 x 10 6 particles for Wo = 8). In the 
colour on-line version of this diagram, simulations using the SPH prescriptions for collision outcome are indicated by red points. Among 
the three models with Wo = 3, N * = 10 8 and Rnb = 0.2 pc, one was computed with the sticky-sphere approximation and experienced 
runaway; the other two, making use of SPH prescriptions, missed the runaway phase. 

This is essentially the same plane as Fig. 1 of Paper I except that i?NB is used instead of R h and that the plotting domain is shifted to 
higher masses and smaller sizes to match the conditions that can be treated by the MC approach and may lead to runaway evolution. 
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computed with the same code and parameters but using dif¬ 
ferent realisations of the cluster and random sequences. In 
three runs, a star grew to 400 — 500 M@ before it left the 
MS; in the fourth, only ~ 300 Mq was reached. To deter¬ 
mine the time at which the runway started, t ra , we look for 
an increase of mass of a factor 3 or more in a star more mas¬ 
sive than 0.9M max (the maximum mass in the IMF) within 
the last tenth of the elapsed time. This definition may seem 
contrived but is required to capture the onset of runaway 
as it would be naturally identified by inspection of collision 
history diagrams such as Figs ll2l or ll3l In practice t ra is very 
close to the core collapse time one would identify by looking 
at the evolution of the Lagrange radii. 

We followed the merger sequence to 1000 Mq at least 
in most cases but, except for a few exceptions, did not try 
to carry on the simulation until the growth was terminated 
by evolution off the MS. With no or little collisional mass 
loss, each merger brings new hydrogen to the VMS and al¬ 
lows it to survive until the next merger, see Figs [H and 
m It is likely that, in real cluster, the growth will saturate 
through some process not included in ME(SSY)**2. One 
possibility is the depletion of the “loss cone”, i.e., the dis¬ 
appearance of the stars populating the orbits that intersect 
the small central volume in which the VMS move around, 
see Section ^. 5. 31 Another limiting mechanism could be that 
the VMS cannot radiate the energy dumped into it by col¬ 
lisions fast enough to keep its relatively small MS size, and 
instead swells and becomes “tr ansparent” to impactors, as 
suggested, for other reasons, bv lGolgatel lll967f) . On the MS, 
a VMS is not only dense enough to stop impactors but also 
resilient to disruption despite its radiation-dominated inte¬ 
rior. Indeed, the binding energy per unit mass actually in¬ 
creases from ~ 3 to ~ 4.5 GM , ^ /R q for masses ranging from 
100 to 10 4 Mq I Bond et al .1984). In any case, we are here 
mostly concerned with detecting the onset of the runaway 
growth and leave the difficult prediction of its final mass for 
future studies. 

Our work is limited to clusters with a relatively large 
number of particles: N * > 10 s , 3 x 10 s , for Wo = 3, 8, 
respectively, because the MC code can only yield robust re¬ 
sults when there are always at least a few dozens particles 
in the cluster core; otherwise, one cannot make a robust es¬ 
timate of the central density, which is required to simulate 
relaxation and compute the collision probabilities. 

An examination of Fig. Q shows that, in general, the 
simple expectation derived from the results of GFR04 that 
VMS formation will occur when and only when the central 
relaxation time is shorter than 20 Myr is borne out by the 
present simulations. However, for relatively small clusters 
(IV* < 3 x 10 5 for Wo = 3; IV* < 10 6 for Wo = 8), a 
shorter central relaxation time appears to be required. This 
is possibly a consequence of the dearth of massive stars in 
the central regions: with our standard IMF, the number of 
stars with M* > 50 Mg in the core of a IV* = 3 x 10 5 , Wo = 3 
or a IV* = 10 6 , Wo = 8 cluster is ~ 30 or ~ 20, respectively. 
Although, over the time required for core collapse ~ 100 Mq 
stars may come to the centre from distances much larger 
than the core radius, the dynamics of core collapse is driven 
by a small number of particles and the applicability of the 
standard treatment of relaxation (and of the notion of 2- 
body relaxation itself) becomes questionable. 


For large clusters, in contrast, we see that runaway VMS 
formation may happen at sizes for which t IC ( 0) > 20 Myr 
because collisions occur from the beginning and accelerate 
core-collapse. Fig. [5] is a graphical attempt at quantifying 
this effect. On this diagram, we plot, for all simulations 
with our standard IMF which lead to VMS formation, £ ra , 
the time of runaway occu rrence. Using equation [8-123] of 
iBinnev fc Tremaine! dl987ll , we define the initial central col¬ 
lision time for a star of 120 Mq as 


Icoll ,C (0) 


tcoll I 


Mi 


where 


A ( N, \ 

- 3/2 , 

'r,\- 7/2 



Up c J 

1 + B 

( Rh > 

V 1 p c I 

K&r 1 


4.13 x 10 12 yr / (M*) \ " 1/2 / Ri + R 2 \ ~ 2 

c n c a ^ m 0 y v r q ) ’ 

22.2 /Ri+R 2 y 1 /Mi + M 2 \ 

c 2 1 Rq J V (M*) )' 


n{0)Rl f 0.385 if Wo = 3, 
IV* ~ 19.10 if Wo = 8, 


CJfj — (ji 


Rh 


GM C i 


0.474 if Wo = 3, 
0.495 if Wo = 8. 


(3) 


We take Mi = 120 Mq, M 2 = (M*) ~ 0.69 Mq and 
Ri ,2 = 77* (Mi, 2 ) according to our MS M-R relation. When 
collisions are treated as perfect mergers (“sticky-sphere” ap¬ 
proximation) , a reasonable fit to our results for the runaway 
time is 


U a _1 = (0.12 trc(0)) -1 + (250 t C oii,c(0)) -1 . (4) 

Hence t ra is approximately 0.12f rc (0), close to the core¬ 
collapse time found in GFR04, 0.15t rc (0), for all clus¬ 
ters with a relatively long initial central collision time. As 
tcoii,c(0) decreases, so does f ra , first mildly and then steeply 
when t C oii,c(0) < 10~ 3 f rc (0). For very short t C oii,c(0), one 
has f ra ~ 250 tcoii,c(0). In this regime, the core collapse ap¬ 
pears to be driven by collisions rather than relaxation. This 
is made clear by comparing the results of two simulations 
of a cluster with Wo = 3, IV* = 10 s (IV P = 3 x 10 5 ) and 
Rnb = 0.05 pc. The first run (K3-52) was realised with the 
usual physics, including 2-body relaxation and collisions but, 
for the second, relaxation was switched off. The evolution of 
the Lagrange radii is plotted for both runs in Fig.[(j] With re¬ 
laxation, core collapse occurs at t cc — 7.4 x 10 -4 Tfp- When 
only collisions drive the evolution (run K3-53), it is slightly 
slower, with f cc — 8.9 x 10™ 4 Tfp, but is otherwise similar. 

It is somewhat surprising that the transition to 
collision-dominated core collapse occurs at such a short 
value of tcoii,c(0) and that t r a/tcoii,c(0) is so long in this 
regime. To some extent, this is due to our definition of 
fcoii,c(0), based on the most massive stars in the IMF (Mi = 
120 Mq), a choice which may not accurately capture the na¬ 
ture of the cluster evolution in the collisional regime. Indeed 
stars much less massive also experience collisions from the 
beginning. Using Mi = 1 Mq, R \ = 1 Mq to define t C oii,c(0), 
one finds t ra — 2t co ii, c (0) when t co n, c (0) <C t rc (0). In any 
case this relation for f ra only holds for sticky-sphere col¬ 
lisions. When the SPH-inspired prescription is used (with 
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mass loss and merger requiring d m i n < Ri + R 2 ) and for 
high velocity dispersions, collisions are less effective at driv¬ 
ing core collapse and f ra is consequently longer for a given, 
small i C oii,c(0) value, as is apparent on Fig. [5] 

Although it is clearly only approximate, it is tempting 
to use equation 0 to predict which conditions will lead to 
collisional runaway by assuming this happens if, and only if 
t ra < 3 Myr. Therefore we have also plotted in Fig01ines in¬ 
dicating this “corrected” runaway condition. One sees that, 
when the velocity dispersion is sufficiently smaller than V*, 
this condition reduces to the condition on the central re¬ 
laxation time alone 1 while at high velocities, one predicts 
runaway for larger cluster sizes than suggested by t TC . How¬ 
ever, the simulation results show the runaway domain for 
very massive clusters to reach still slightly larger f?NB val¬ 
ues. Some of this may be due to dispersion in the f ra results 
(see below) but, by inspecting the run K3-55, another effect 
is discovered. In this case, the VMS star survives to at least 
t = 4 Myr, despite its MS lifetime being shorter than 3 Myr. 
This is due to the combination of (1) the prescription used 
to set the effective age of merger products by following their 
central He content (“minimal rejuvenation”, see Paper I), 
(2) the relation we use for the value of the core He at the 
end of the MS as a function of the stellar mass (Fig. 3b of 
Paper I) and, (3) for this particular simulation, the assump¬ 
tion of perfect merger. Because the MS lifetime is nearly 
independent of the mass for M* ^ 100 Mg but more mas¬ 
sive stars transform a larger fraction of their mass into He 
on the MS, each merger amounts to an effective rejuvena¬ 
tion if no mass-loss is allowed. Clearly collisional mass loss 
may change the picture in this very special regime but we 
think the most important shortcoming of the approach is to 
assume the VMS to be on the MS. As we will see below, 
in most cases, the average time between collisions is much 
shorter than the thermal timescale of a MS VMS so that it 
structure is likely strongly affected by collisions and our SPH 
mass-loss prescriptions would not be appropriate anyway. 

The extension of the runaway domain to larger sys¬ 
tems opens the possibility that proto-galactic nuclei with¬ 
out MBH may be subject to runaway formation of a cen¬ 
tral massive object, despite relatively long relaxation times. 
Present-day nuclei have a size of a few pc. However, using 
equation 0 as a guide, a compact J?nb = 1 pc nucleus with 
Wo = 8 needs to be at least as massive as ~ 3 x 10 9 Mq to 
experience collision-driven collapse in less than 3 Myr. For 
_R nb = 3 pc, the minimum mass is ~ 3 x 10 10 Mg! So it seems 
that this process may operate in galactic nuclei only if they 
are born very compact or concentrated. Indeed for Wo = 10, 
one would expect all clusters more compact than ~ 2.5 pc 
to complete core collapse within 3 Myr, independent of their 
mass. 

Our results show a relatively large dispersion in the 
runaway time when evolution is dominated by relaxation, 
in which case t ra should be very close to the core-collapse 
time found in GFR04, i.e., tcc| rlx — 0.15f rc (0) (with little 
variation from one simulation to another). The present re¬ 
sults indicate a slightly faster evolution, t cc | r i x — 0.12t rc (0), 
which is very likely the consequence of the use of a different 


1 This had to be expected because the ratio of relaxation time to 
collision time is t r i x /t co n ~ (In A) - 1 (l+0)0 —2 with 9 — (V*/cr v ) 2 . 


Time [Myr] 

0.01 0.1 1 10 10 2 



Figure 3. Evolution of Lagrange radii for model K8-27. For this 
cluster, the core collapse is slightly too long to allow the runaway 
process to set in. The cluster undergoes two successive core col¬ 
lapses. The first one driven by the massive MS stars, the second 
by the stellar BHs (with a mass of 7 Mq) resulting from their 
evolution off the MS. For clarity, the curves have been smoothed 
by using a sliding average procedure with a (truncated) Gaussian 
kernel over 10 adjacent data points. 

See also Fig.Q] 

code 2 . The larger intrinsic dispersion is also a property of 
the present code, linked to the use of local time steps; they 
impose a scheme in which the selection of the next particle 
pair to evolve is a random process. 

Following this overview of our results for the complete 
set of standard simulations, we proceed, in the following 
subsections to a more detailed description of the runaway 
process based on a few illustrative cases. Presentation of 
non-standard cases, for which we varied the IMF, M-R re¬ 
lation, collision prescription, or the treatment of the VMS, 
will be done in Section HU 

2.3 Cluster structure evolution 

2.3.1 Missing the collisional runaway 

Figs EH and 2] illustrate what occurs when the core-collapse 
time is longer than MS lifetime of massive stars. For the 
sake of clarity, we present a case in which the segregation- 
induced core collapse would have taken just slightly longer 
than t* = 3 Myr. The evolution of the Lagrange radii, plot¬ 
ted in Fig. [I] indeed indicates a first core contraction which 
stopped and reversed at t*. A second collapse occurs much 
later, at t ~ 500 Myr. What happens is made clear by 

2 See [Freitag_^^en3 I 2 OO ll) for comparisons of the tcc| r i x val¬ 
ues found with different codes, in the case of the collapse of a 
single-mass model. Models with a broad mass spectrum are more 
affected by particle noise than single-mass clusters (at a given 
N p ) and display more intrinsic and code-to-code variation. 
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Time [Myr] 


Figure 4. Evolution of the various stellar types for model K8-27 
(same as Fig. E3 For this cluster, the core collapse is slightly too 
long to allow the runaway process to set in. The stellar types 
accounted for are: MS stars (MS), white dwarfs (WD), neutron 
stars (NS) and stellar BHs (BH). For simplicity, all WDs have 
0.6 Mq, all NS 1.4 M© and all BHs 7M©. Top panel: number of 
compact remnants as function of time. Bottom panel: evolution 
of Lagrange radii. The simulation was stopped as the stellar BHs 
underwent core collapse on their own. 


Fig. a On the top panel, we have plotted the evolution of 
the number of stellar remnants; the bottom panel is the 
evolution of the Lagrange radii for each stellar type (MS: 
main-sequence stars, WD: white dwarfs [0.6 Mq], NS: neu¬ 
tron stars [1.4 Mq], BH: stellar black holes [7Mq]). The first 
collapse stops when the massive stars turn into BHs, causing 
strong mass-loss from the central region which re-expands 
as its binding energy decreases. The BHs are born strongly 
segregated because their massive progenitors had concen¬ 
trated at the centre. The BH distribution first re-expands 
but eventually re-collapses as a result of mass segregation. 
From i ~ 40 Myr onwards, they are indeed the most mas¬ 
sive objects in the cluster. The subsequent evolution of this 
population of BHs will be driven by their binaries, either 
present before the second collapse, or dynamically formed 
through 3-body interactions during it. The MC code used 
here cannot treat binaries so we leave the study of the evo¬ 
lution of t he central dens e BH sub-cluster out of the present 
work fsee lO ’Learv et HM). 


2.3.2 Core collapse with runaway 

The evolution of Lagrange radii during the core collapse and 
runaway VMS growth are depicted, for four Wo = 3 cases, 


Time [Myr] 



Figure 5. Evolution of Lagrange radii for various clusters with 
Wo = 3. From the least to the most collisional (initially): with 
solid lines (black in colour version), K3-33; with long dashes (blue) 
K3-37; with short dashes (magenta) K3-64; with dotted lines (red) 
K3-61. The top horizontal and right vertical axis indicate time 
and radius in physical units for model K3-33. The evolution of 
this cluster is driven by 2-body relaxation until late into core¬ 
collapse, at which point the collisions kick in. The core collapse 
time is thus very similar to the value one would obtain for point- 
mass dynamics. 


Time [Myr] 

0 0.05 0.1 



Figure 6. Evolution of Lagrange radii for model K3-52. This 
cluster is strongly collisional from the beginning. The core col¬ 
lapse is driven mostly by mergers and occurs on a timescale much 
shorter (when measured in FP units). The dashed lines show the 
evolution of the same system without relaxation (K3-53). Core 
collapse still occurs, although on a slightly longer timescale. 
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Figure 7. Evolution of the density profile for all stars with 
masses between 0.2 and 120 M® (i.e., excluding the runaway ob¬ 
ject) in simulation K8-29. For each profile in the sequence, we 
indicate the time and mass of the most massive star in the clus¬ 
ter at this time. The FP time unit is Tpp ~ 4.72 x 10 9 yr. A 
central sub-cluster with a cuspy density profile forms. Its profile 
is compatible w i th a p ower-law of exponent —1.75 but it is not a 
iBahcall &; WolJ Jl97ol cusp (see Figs 181 FTT)1 and text). 



Figure 8. Similar to Fig. Qbut we only plot the density of stars 
with masses between 3 and 120 Mq. This figure confirms that the 
most massive stars are responsible for the central detached cusp 
seen in Fig. [3 However, it seems that the innermost part of the 
density profile is more concentrated than a p oc R ~ 175 cusp (see 
Fig.[HJ. 


Figure 9. Density profiles for various stellar masses in the same 
simulation as in Figs Q and [s] To investigate the structure of the 
cluster during the runaway phase, we have averaged the density 
profiles for 18 snapshots within the given time interval. The var¬ 
ious curves correspond to different ranges in stellar mass. It is 
clear that only stars more massive than 10 Mq (with a number 
fraction of only ~ 5 x 10 —3 ) contribute to the formation of the 
cusp and that its profile is steeper than p oc R ~ 1 - 75 . The line with 
slope -2 is plotted as a guide. Note that the radius and density 
scales and ranges are different from those of Figs □ and 0 


in Fig. 0 The first simulation is for a cluster which not is 
collisional initially, i.e., from Eq. Q] with i C oii,c(0) —as de¬ 
fined in Eq.|3for 120Mq stars— larger than ~ 10~ 3 f r c(0): 
£coii,c(0) = 3.0 x 10 6 yr = 0.056 t rc (0). For this particular 
cluster, core collapse occurs in approximately 2.5 Myr, thus 
beating stellar evolution. Collisions are very rare until the 
runaway starts in deep core collapse, see Fig. m Therefore, 
the evolution to core collapse is very similar to what was 
obtained for the pure relaxation case in Paper I. Once the 
runaway has started, the central regions re-expand as a re¬ 
sult of the heating effect of the accretion of deeply bound 
stars by the central object, as in a cluster with a ce ntra l BH 
destroying or ejecting stars that come too close iShapirol 


1977]: (Duncan & Shanirol 

19821: lAmaro-Seoane et alJ 120041: 

Baumsardt et al. 12004a b 

). We cannot follow this “post- 


collapse” evolution very far because our code slows down 
considerably when reaching the runaway phase, as the cen¬ 
tral relaxation and collision times plummet. Furthermore, 
as will be discussed in Section ^. 5. 31 it is likely that our way 
of computing collision probabilities applies correctly to the 
VMS only at the beginning of its growth. 

The other simulations plotted in Fig. [3 correspond 
to denser clusters in which collisions play an increasingly 
important role in the dynamics and lead to shorter and 
shorter collapse time (when measured in relaxation times). 
The central collision times t C oii,c(0) for a 120 Mq star are 
6.2 x 10 5 yr = 0.014 t rc (0), 8.8 x 10 3 yr = 1.3 x 10" 4 t rc (0) 
and 500yr = 2.9 x 10~ 5 t rc (0). As explained in Section [2.21 
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Figure 10. Similar to Fig. [H] but for the 1-D velocity disper¬ 
sion a v (plotted on a linear scale). As a result of (partial) energy 
equipartition, the most massive stars are considerably cooler than 
other stars in the central regions. The dotted line is the initial ve¬ 
locity dispersion (Wo = 8 King model). The dash-dotted line 
indicates the Keplerian dispersion ctk = y/ GMy ms /(3 R) corre¬ 
sponding to a central object with My ms = 1000 Mg. One sees 
that the sphere of influence of the runaway star (inside which 
<r„ ~ ctk) is not resolved. 



Figure 11. Evolution of the central density (p c , dashed line, left 
scale) and 3-D velocity dispersion (as c , solid line, right scale) for 
model K3-37. Note how the central parts cool down during core 
collapse as a result of partial equipartition between the massive 
stars and the lighter ones. 


in extremely dense clusters, core collapse is driven by colli¬ 
sions themselves, with relaxation playing only a minor role. 
This is illustrated in Fig. |fi] for another Wo = 3 cluster. 
For these runs, we have assumed that all collisions result in 
complete mergers. In contrast to gravitational 2-body en¬ 
counters, such mergers do not redistribute energy between 
stars; they continuously dissipate orbital energy, leading to 
shrinkage of all inner Lagrange radii with no rebound. 

To get a more precise description of the structure of a 
cluster while it undergoes core collapse and runaway VMS 
formation, we look at successive radial profiles of the stellar 
density and velocity dispersion. To have a sufficient resolu¬ 
tion of the central regions, which are the most affected by 
the dynamical evolution, we use our highest resolution simu¬ 
lation, model K8-29, with 9 x 10 6 particles (the largest num¬ 
ber that would fit into the memory of available computers). 
The initial structure is a Wo = 8 King model. Fig. Q shows 
the evolution of the total stellar mass density. The core col¬ 
lapse is characterised by the formation a small central den¬ 
sity peak, growing inside the initial constant-density core. 
In comparison with the core-collapse of single-mass models, 
the evolution shown here is relatively unremarkable and does 
not appear to proceed in a self-similar manner. In particular, 
the contracting cor e does not leave behind a p oc R 1 re¬ 
gion with 7 ~ 2.23 iFreitae fc Bendl 200 ll : iBaumeardt et alJ 
120031 and references therein). The central density profile in 
deep core collapse is better fitted with 7 ~ 1.75. It may be 
tempting to interpret this value of the exponent by th e for¬ 
mation of a Bahcall-Wolf cusp feahcall fc Wolfll97fiil , sim¬ 
ilarly to what 2 -body relaxation creates in a stellar system 
dominated by a central massive object. This would not be 
correct, however, as the following figures make clear. When 
we look at the evolution of the distribution of stars with 
masses in the range between 3 and 120 Mq in Fig. [H] we 
see a much stronger evolution than for the overall distribu¬ 
tion because lighter stars, which form an essentially static 
background, are so much more numerous. Fig. [§] yields a 
detailed picture of the distribution of stars in various mass 
ranges during the late stages of the core collapse (while the 
VMS is quickly growing). Clearly, only the spatial distribu¬ 
tion of stars more massive than 10 Mq significantly changes 
and the central density peak is mostly made of stars more 
massive than 30 Mq. These objects do not form a 7 ~ 1.75 
cusp but rather a distribution of varying slope. 

In Fig. EH are plotted the velocity dispersion profiles 
for the same ranges of stellar masses as in Fig. EJ Most im¬ 
portantly, this figure shows that, inside a region roughly 
corresponding to the initial core (R c ore(0) = 0.12 pc), relax¬ 
ation has produced partial energy equipartition with stars 
above 10 Mq having cooled to lower velocities by heating 
up lighter stars. This is particularly dramatic for the stars 
in the 30 — 120 Mq bin, which exhibit an inverse tempera¬ 
ture gradient. Around R = 0.2 the velocity dispersion for 
these massive stars has in fact increased above the initial 
value, because the stars with higher orbital velocity are 
less affected by dynamical friction and have stayed in place 
while others were sinking. In this figure, we also plot the 
Keplerian velocity for a Mvms = 1000 Mq central object, 
(jk = y / GMvms / (3 R). Clearly, we cannot resolve the radius 
of influence (the central region dominated by the central ob¬ 
ject’s gravity) for any range of stellar mass. Only within the 
radius of influence would one expect a Bahcall-Wolf cusp to 
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Figure 12. Collisional history for run K3-33. We represent the 
evolution of the five particles that have experienced the largest 
number of collisions. The top panel represents the radius (distance 
from the cluster centre) where the collision occurred. The bottom 
panel shows the evolution of the mass. A circled dot indicates 
that the star has merged with a larger one (usually the runaway 
star). Note that there are few collisions until t — 2.5 Myr, i.e., the 
moment of core collapse, at which time one star starts growing 
very quickly. 

develop. The evolution of the central density and velocity 
dispersion (for another simulation) is depicted in Fig. 1111 
We see how the central density steadily increases while (in 
contrast with the single-mass core collapse) the velocity dis¬ 
persion drops, a result of energy equipartition. Both effects 
contribute to making collisions more and more likely. Fur¬ 
thermore, the decrease in velocity dispersion ensures that 
disruptive collisions will be unlikely unless the central ve¬ 
locity dispersion is initially very high. 

2.4 Collisional runaway sequences 

We first describe results obtained in the sticky-sphere ap¬ 
proximation for collisions. We will then present the differ¬ 
ences introduced by the more realistic collision treatment 
based on SPH results. 

A typical runaway history for a cluster which is not 
initially collisional is depicted in Fig. ini where we follow 
five stars participating in collisions, including the runaway 
merger product. The parameters of the cluster are such that 
it undergoes core collapse in ~ 2.5 Myr. Only a few collisions 
happen before deep core collapse, at which time one star de¬ 
taches itself from the rest of the population by growing in 
a rapid succession of mergers. Two characteristics of this 



Figure 13. Collisional history for run K3-52. Collisions occur 
from the beginning of the evolution and many stars grow in an 
orderly way until one detaches itself from the others and starts 
growing in a runaway fashion. See text for explanations. 

runaway growth are its steepness and the fact that no other 
star experiences significant mass increase. The initial cen¬ 
tral velocity dispersion in this cluster is a o ~ 52kms _1 , a 
value for which we expect the sticky-sphere approximation 
to be accurate. Fig. cna is an example of runaway in a clus¬ 
ter which is initially strongly collisional. In this case, many 
stars experience numerous collisions and start growing from 
the beginning of the evolution. Nevertheless, only one of 
these objects experiences a very fast growth at the moment 
of (collision-driven) core collapse. Here cro ~ 1200 km s -1 
but the escape velocity from a 1000 Mq VMS on the MS 
is Kbc — 3000kms _1 , so gravitational focusing is still ef¬ 
fective, producing the steep increase of cross section with 
My ms required for runaway. Because relative velocities are 
still significantly lower than 14 sc , we do not expect strong 
collisional mass loss but the critical d m i n for merger becomes 
significantly smaller than Ri + R 2 , making mergers less fre¬ 
quent than assumed here. 

To study which stars contribute to the build-up of the 
VMS, we draw merger trees. Four examples are presented in 
Figs 1141 and 1151 The first is for two clusters with Wo = 3 
and IV* = N p = 3 x 10 5 , the second for Wo = 3 and 
IV * = N p = 3 x 10 6 . In each figure, panel (b) shows the 
case of a smaller, initially more collisional cluster. In all 
cases, we follow the VMS growth up to 2000 Mq. The clus¬ 
ter in Fig. im has .Rnb = 0.4 pc, leading to collapse in 
tec = 2.9 Myr, just in time to produce collisional runaway 
before the massive stars would have left the MS. In this 
situation, typical of clusters which first become collisional 
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(a) 


(b) 


Figure 14. Merger trees for simulations of 
clusters with Wo = 3, TV* = N p = 3 x 10 5 . 
Panel (a): -Rnb — 0.4pc (model K3-23). 
Panel (b): Rnb — 0.03 pc (model K3-9). 
We follow the growth of the runaway star 
to ~ 2000 Mq. The horizontal axis indi¬ 
cates the mass of stars taking part in the 
collisions. The sequence of mergers is rep¬ 
resented in the vertical direction from top 
to bottom. The right “trunk” is the grow¬ 
ing VMS. The radius of the disks is propor¬ 
tional to that of the corresponding star. In 
case (a), the cluster is not initially colli- 
sional (the collision time for a 120 M© star 
is much larger than its MS life-time) and 
most collisions occurs in deep collapse and 
feature stars of mass 60—120 M© that have 
segregated to the centre. Most of them are 
not themselves merger products. In case 
(b), the cluster is initially collisional and, 
although the runaway also occurs in core 
collapse, most stars contributing to it have 
experienced earlier collisions. 




Figure 15. Merger trees for the simu¬ 
lation of a cluster with Wo = 3, TV* = 
Np = 3 x 10 6 . Panel (a): 1 ?nb — 0.2pc 
(model K3-37). Panel (b): Rnb = 0.1 pc 
(model K3-36). We follow the growth of 
the runaway star to ~ 2000 M©. Both clus¬ 
ters have initial central collision times for 
120 M© stars shorter than 3 Myr, hence the 
relatively complex merger trees. The more 
compact cluster, case (b), is more colli¬ 
sional which explains a stronger contribu¬ 
tion of low-mass stars. 
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Figure 16. Evolution of mass and various timescales during the runaway for the same simulations as in Fig. rm We plot the time between 
successive collisions, an estimate of the Kelvin-Helmholtz timescale ^kh °f the collision product (assuming normal MS M-R relation 
and luminosity), the time left until exhaustion of hydrogen at the centre (from our simple “minimal rejuvenation” prescription), as well 
as the mass of the star (right scale). No use is made of £kh during the simulations. We plot it here to show than in the vast majority of 
cases, the interval between collisions is (much) shorter than £kh so that the star would probably be out of thermal equilibrium, with a 
swollen structure. 
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Figure IT. Mass and timescales evolution during the growth of the runaway star for the same simulations as in Fig. m See caption of 
Fig. 1101 for explanations. 


in deep collapse, the merger tree has a very simple struc¬ 
ture, with most branches lacking sub-structure: the stars 
that merge with the VMS (the trunk) are not themselves 
collision products. Furthermore, a large majority of these 
impactors have masses around the tip of the IMF, between 
60 and 120 Mq . These stars add up to 92 % of the VMS 
mass. The distribution of the number of contributing stars 
(the tips of the branches, i.e., stars that are not themselves 


collision products) is bimodal with one peak between 0.2 
and 1 M 0 and another covering the 40 — 120 M 0 ; it reflects 
the MF in the central regions (GFR04). This is to be con¬ 
trasted with Fig. m>, obtained for a much more collisional 
cluster with Rnb = 0.03 pc, which enters the runaway phase 
at t ~ 0.064 Myr. Here most stars merging with the VMS 
were themselves produced by previous mergers and the con¬ 
tribution of lighter stars is larger, although the 60 — 120 M 0 
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range still contributes 78 % of the VMS mass. The cause of 
the difference is that, in the more compact cluster, collisions 
occur from the beginning, when the cluster is not segregated 
yet; in the larger cluster, virtually all collisions take place in 
the high-density core formed by the concentration of high- 
mass stars. The two clusters in Fig. 1151 with 3 x 10 6 stars 
and sizes of 7 ?nb = 0.2 pc and 0.1 pc are both relatively colli- 
sional, as required for this Wo value and this number of stars 
to experience core collapse in less than 3Myr (see Fig. 0: 
thus, they exhibit complex merger trees. Again, in the more 
compact cluster, low-mass stars play a larger but not pre¬ 
dominant role, with some 90 % of the VMS mass originating 
in stars more massive than 60 M@. Merger trees for Wo = 8 
clusters exhibit the same general characteristics and variety. 
In particular, the same kind of bimodal distribution in the 
masses of the contributing stars can be observed, with more 
stars in the 0.2 — 1 Mq range for clusters with shorter initial 

1'col [ /trlx • 

It is instructive to examine the various timescales in¬ 
volved in the runaway process. This is done in Figs Em 
and EH for the same four simulations as in Figs rm and 
EH In these diagrams, we follow the VMS, merger after 
merger, making use of two ordinate axis, one to indicate 
its mass, the other to monitor important timescales. Those 
are the time to the next collision Af co n, the time left on 
the MS AtMS (if the star were left to evolve without fur¬ 
ther collision), and the MS thermal relaxation time, i.e., the 
Kelvin-Helinholtz (K-H) time. The latter is approximated 
here by £kh = GM*(J?*L *) _1 where L* is the luminosity 
of the star (c l ose to the Eddington limit for VMS; see, e.g., 
iBaraffe et alJl2f)0ll : ISchaererlteOOSf) . From 100 to 10 4 Mq, 
£kh increases slowly from ~ 10 4 to 6 x 10 4 yr. Here we do 
not truncate at My ms = 2000 Mq , but show instead the full 
extent of the simulations. However, in all but the first case 
(Fig. ESK)> we stopped the simulation before the VMS left 
the MS, when the average interval between mergers was still 
much shorter than the remaining MS lifetime. 

Because it corresponds to a core-collapse time just 
slightly shorter than 3Myr, the first simulation of the four 
is one of the exceptional few for which we could afford to 
integrate until stellar evolution caught up with collisions, 
as indicated by the convergence of the At co ii and A£ms 
curves. In most other cases, A£ms still exceeded A£ms by 
orders of magnitude —suggesting many more collisions to 
come - when we decided to terminate the simulation. From 
the merger trees and the curve showing My ms as function 
of the collision number, there is only little sign of a decrease 
with time, as massive stars that have segregated to the cen¬ 
tre are progressively depleted. Such trend can be seen for 
the runs with 3 x 10 5 stars (in particular if one draws the 
complete merger tree past My ms = 2000 Mq ), but is not 
perceptible for a cluster of 3 x 10® stars whose central parts 
constitute a larger reservoir of massive stars. Indeed among 
3 x 10 s stars forming a Wo = 3 cluster, there is only about 
1600 Mq of stars more massive than 60 Mq in the core. 

An essential result is that, in the vast majority of cases 
the K-H time is (much) longer that the average time between 
successive collisions. This means that the VMS is unlikely 
to relax to thermal equilibrium, i.e., to the MS after it has 
merged with an impactor and before the next collision. In 
principle, then, it may be kept in a swollen state by re¬ 
peated collisions and our assumptions of MS M-R relation 



Figure 18. Collision parameters for a cluster with relatively low 
velocity dispersion (same simulation as in Figs ll5h , and 117b - ). For 
each collision (dot), we indicate the “impact parameter”, in units 
of the stellar half-mass radii, and the relative velocity at large 
separation, in units of see Eq. 13 of Paper I). Circled dots 

indicate mergers involving a star with a mass larger than 120 Mq 
( i.e., the runaway object). The lines indicate the condition for 
merger (Eq. 13 of Paper I) for [Mi,M 2 ]/Mq = [1,1] (solid line), 
[1,10], (dots), [0.1,1], (short dashes), [0.1,60], (long dashes) and 
[60,60], (dot - dash). Only below these lines should the collisions 
result in mergers. 


and continuous nuclear burning (to set A£ms) are probably 
incorrect during runaway growth. This suggests a very dif¬ 
ferent picture concerning the evolution of the VMS in the 
runaway phase and the termination of the growth. The VMS 
may grow larger and larger until it is so diffuse that it lets 
the smaller stars fly through it without the ability to stop 
them, h ence bringing back the “transpa rency saturation” 
jColeatelll967f : iLiehtman fc Shapirdll978i l. Another impor¬ 
tant effect is that the central nuclear reactions are very likely 
switched off in a merg er product which expands as a result 
of shock heating le.g.. iLombardi et al.ll200^ . 1200311 . so the 
stellar evolution of the VMS may be suspended. Of course 
these effects will persist only as long as A£ co ii < £kh- When 
the collisions eventually become rarer —from exhaustion of 
stars to collide with or because stellar evolution from nor¬ 
mal massive stars drives core re-expansion at t ~ 3 Myr— 
the VMS returns to the MS and to nuclear burning. Given 
the complexity of interplay between repeated collisions and 
VMS structure (itself determining the probability and out¬ 
come of further collisions), it seems impossible to make any 
simple prediction as to what the final VMS mass would be 
when it eventually leaves the MS. In future works, this situ¬ 
ation may be studied by combining a stellar dynamics code 
with on-the-fly SPH sim ulation to fo ll ow the VMS stru cture 
collision after collision feaumgardt fc Nakasatci 12005 ). For 
the time being, these aspects have to be considered as one 
central source of uncertainty regarding the outcome of the 
runaway phase. 
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Figure 20. Role of mass-loss in the runaway process. For the three simulations presented here, we used the same initial model with 
Wo = 3, TV* = 3 x 10 5 and Rnb — 0.3pc (models K3-13, K3-14 and K3-15). All collisions are assumed to result in merger with a 
constant fractional mass loss 5m = 5M/M , i.e., M me rger = (1 — 5m){M\ + M 2 ) where Mi, M 2 are the masses of the colliding stars. This 
prescription is artificial in that a collision with a light star will cause more damage than one with a more massive star. For case (a), with 
no mass loss, the simulation was interrupted before the runaway object evolved off the MS because mergers constantly rejuvenate the 
star. With some mass loss, stellar evolution can catch up with collisional growth and terminate it. For simplicity, in these cases, we did 
not assume an IMBH was formed at the end of the VMS MS, but a 7 Mq BH instead, which was not allowed to collide with other stars. 
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Figure 19. Same as Fig. ll8l but for a cluster with high velocity 
dispersion: model K3-55. For this case, a large fraction of collisions 
would actually be fly-bys and the assumption of perfect mergers 
is questionable. 


Figure 21. Cumulative distribution of fractional mass loss in 
collisions for four cluster simulations using the SPH prescription. 
In all four cases, runaway VMS formation happened. The average 
fractional mass losses for these runs are 0.2, 0.4, 0.5 and 4.9%. 


Putting aside the possible complications due to the non- 
MS structure of the VMS, we can look at the changes in the 
runaway process introduced by a more realistic treatment 
of collisions between MS stars, as allowed by the SPH re¬ 
sults. Figs EU and 1191 show when and why the sticky sphere 
approximation may become questionable. These diagrams 


show the impact parameter and relative velocity for all col¬ 
lisions that happened in two cluster simulations, one with 
ao ~ llOkms -1 CFig. 1181 . the other with ao — 635 km s 1 
(Fig. 1 1 91 . For lower velocity dispersion, a large majority 
of encounters happened with parameters leading to merger 
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according to Eq. 13 of Paper I. In particular, all colli¬ 
sions but one in which the runaway object took part were 
in this regime. They had relative velocities smaller than 
0.1 V t and would have produced less than 1 — 3 % mass loss 
llFreitae fc BendEoOfih . Consequently, treating all collisions 
as perfect mergers is certainly a very good approximation 
for system with comparable or lower velocity dispersion, in¬ 
cluding any reasonable globular cluster. In the high-velocity 
situation, a significant fraction of the collisions (featuring 
the runaway object or not) should have left the stars un¬ 
bound rather than produced a single merged object, sug¬ 
gesting that using the SPH prescription to set the impact 
parameter for merger may affect the results. On the other 
hand, most collisions occurred at less than V* and would 
have produced no more than a few % mass loss. 

To assess the impact of mass loss on the runaway mech¬ 
anism, we have performed a few runs for which we assumed 
that all collisions result in a merger with some constant frac¬ 
tional mass loss 5m, varying the value of 5m- We show run¬ 
away growth sequences for 5m = 0, 0.05 and 0.10 in Fig. 03 
One sees that 5m ^ 0.10 is required to prevent the growth 
of a star more massive than a few 100 M@ . As typical im- 
pactors have a mass around Mi mp ~ 100 Mg , it is obvious 
why a VMS of 1000 Mg can only be grown through a se¬ 
quence of merger if 5m < 0.1. The fact that we assume 
all mass lost to originate from the hydrogen envelopes also 
shortens the remaining MS lifetime of the collision product 
and prevents reaching My ms — M^/Sm- 

SPH simulations show mass loss of at most a few per¬ 
cent for < V t , so reaching 5m ~ 0.1 (on average) would 
require really extreme conditions. The cluster with the high¬ 
est velocity dispersion for which we have performed a sim¬ 
ulation using the SPH prescriptions is model K3-51, with 
(Jo — 1270 kms -1 . It has the same parameters as the sticky- 
sphere run K3-61 and we obtain the runaway VMS growth 
in a qualitatively similar fashion, although about 4 times 
more collisions are required to attain a given My ms- The 
most noticeable difference is that the time required to enter 
the runaway phase t ra is also about 4 times longer than in 
the sticky-sphere approximation (see left-most open triangle 
in Fig. 0 . This is obviously a consequence of the collisions 
being about 4 times less efficient at driving collapse in this 
collision-dominated cluster. In cases with lower velocity dis¬ 
persion, using the SPH prescription make very little change 
to the conditions required to achieve fast collapse and colli- 
sional runaway or to the runaway itself. 

In Fig. 0] we plot the cumulative fractional mass loss 
for all collisions that took place in four simulations for clus¬ 
ters with initial central velocity dispersions ranging from 
(Jo = 28 to 1420kms _1 . We see that all collisions yield 
5m <10% for (Jo < lOOkms -1 with an average (5m) < 1 %. 
The cluster with the highest uo experiences a few completely 
disruptive collisions but only involving stars less massive 
than 2 Mg. Even in this extreme case, (5m) is only of or¬ 
der 5%. For collisions involving a star more massive than 
120 Mg, the average mass loss is as low as 0.3 %. 

The only case for which we find a qualitatively differ¬ 
ent evolution when using SPH-generated collision recipes 
instead of pure mergers, is for Wo = 3, V* = 10 8 and 
Rnb = 0.2 pc. In K3-55, assuming pure mergers, we found 
runway, thanks to collisional rejuvenation that kept the 
VMS on the MS for more than 4 Myr. With collisional mass 



M * in M s 


Figure 22. Merger tree for the runaway growth of a VMS in 
model K3-61, a cluster with steep Kroupa IMF extending from 
0.01 to 120 Mg (see text). 

loss and a large fraction of collisions not resulting in merg¬ 
ers, the (orderly) growth of stars stopped at i ~ 3.5 Myr, 
before any runaway star detached from the mass spectrum 
(runs K3-65 and K3-56). 

Because a realistic treatment of collisions has so little 
impact on most results pertaining to the cluster evolution 
and runaway VMS growth, it is no surprise that we find them 
essentially unchanged if, instead of our SPH-based prescrip- 
tion, we im plement the simple fitting formulae published by 
iRauchl (; 1999fl . which were derived from a smaller set of SPH 
simulations. 

2.5 Non-standard simulations 

2.5.1 Other IMFs 

GFR04 established that, for any realistic IMF and all clus¬ 
ter models we tried for which t rc (0) > 0, the segregation- 
driven core collapse occurs at f C c| rlx — 0.15t r c(0). By “re¬ 
alistic” IMF we mean one that is sufficiently broad, with 
[i = Mt.m^/jM t) > 50, and not too steep. This includes the 
“Kroupa” IMF dKroupa et a 1HD, which contains signif¬ 
icantly fewer high-mass stars than the Salpeter IMF (with 
dW/dM* oc MV 2 ' 7 for M* > lMg). Although this is ac¬ 
tually an effective IMF for stars forming the Galactic held, 
with c lusters probably having a Salpeter-like IMF beyon d 
1 Mg dKroupa fc Weidnerll2003t fWeidner fc Krouoall2004ll . 
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we used it in order to investigate the possible effects of a 
smaller fraction of massive stars. We find that the runaway 
growth of a VMS happens essentially in the same way as 
with the Salpeter IMF. The only noticeable, and not unex¬ 
pected, difference is that because stars around 100 Mq are 
so much rarer, the average impactor mass decreases much 
more quickly during the VMS growth, as can be seen in the 
merger tree shown in Fig. 1221 

We have also simulated clusters with a reduced stel¬ 
lar mass range. The single-mass case lacks realism and was 
only considered in Paper I as a test-case to compare with 
iQuinlan fc Shaoirol (Il99ff) . 

If the Salpeter IMF only extends from 0.2 to 10Mg, 
the average stellar mass is 0.58 Mq and /r = 10/0.58 ~ 17. 
The results of GFR04 suggest t C c| rlx — 0.7 t rc (0) for Wo = 3. 
For such a cluster, we find i ra — 1.2f rc (0) (run K3-20). The 
10 Mq stars leave the MS after t* = 25 Myr; these two nearly 
equally longer timescales combine to give conditions for run¬ 
away in the (W, 7?nb) plane that are very close to those for 
our standard IMF. Hence, the details of the upper end of 
the IMF have little effect on the conditions for the onset of 
runaway collisions. Of course, if the maximum stellar mass 
in the IMF is lower, a larger number of mergers are required 
to grow a VMS but there is about 3 times more mass in 
stars in the range 5 — 10 Mq in a Salpeter IMF truncated 
at 10 Mq than in the range 60 — 120 Mq if the IMF extends 
to 120 Mq. In simulation K3-20, the cluster entered the run¬ 
away regime at f ra ~ i C c| rlx = 20 Myr, allowing a VMS to 
reach 1500 Mq when we stopped the computation. 

2.5.2 Conditions for collisional runaway in MGG-9 

IPortegies Zwart et alJ ([2004) have performed IV-body sim¬ 
ulations of the clusters MGG-9 and MGG-1 1 in M82 using 
~ 130 000 to ~ 600 000 particles. Following iMcGradv et all 
(2003), they assume for MGG-9 an initial mass of 1.8 x 
10 6 Mq and a half-mass radius of 2.6 pc. Unlike the case of 
MGG-11 for which there are indications of a lower cut-off of 
the IMF at ~ 1 Mq (which yields (M») ~ 3Mq), the obser¬ 
vational data for MGG-9 is compatible with a normal cluste r 
Kroupa IMF extending from 0.1 to 100 Mq llKrou pall200Ji : 
dNt/dM „ oc Mf a with a = 1.3 below 0.5 Mq and a = 2.3 
for higher masses ((M,) ~ 0.64 Mq). To get the total cluster 
mass as estimated observationally from its size and velocity 
dispersion, some 2.7 x 10® stars are required. Because this 
numb er is too hig h for present-da y direct IV-body simula¬ 
tions, |Po^gie^Zwm^^dj (2004) adopted a 1 Mq cut-off. 
They considered King models with Wo ranging from 6 to 12; 
none of them led to collisional runaway, in contrast to their 
MGG-11 models. 

In light of GFR04, this is a surprising result. From its 
position in Fig. 1 of Paper I, the MGG-9 model is expected 
to undergo quick core collapse and runaway collisions if its 
initial concentration corresponds to Wo >8.5. We have run 
three simulations of this cluster with N p = V* = 2.7x 10 6 for 
Wo = 8, 9 and 12 (MGG9-K8, MGG9-K9 and MGG9-K12a). The 
model with the lowest concentration does not experience 
runaway but the other two do. We did not try to simulate 
MGG-11. If the average stellar mass is as high as 3Mq, the 
number of stars in that cluster is of order ~ 10 s , too small 
for robust MC treatment with a broad IMF and Wo > 3. 
We note, however that the condition for runaway suggested 



Figure 23. Conditions for collisional runawa^Tommodelsofthe 
clusters MGG-9 and MGG-11 in M82 l MglltTulv_eL_aI] [2111131) - 
This diagram is adapted from Fig. 3 of^^^^^^^^^uTeTaLl 
l2004f) . Their dynamical friction time to f is approximately equal 
to 4.7 x 10“ 3 Tpp for models with (M*) = 3Mq and 1.3 X 
10 -3 Tpp for (M,) = 0.64Mq. c is the initial concentration pa¬ 
rameter; corresponding Wo values are gi ye n on the right of the di- 
agram. Round points are for simulations by IPort ogies Zwnrt et al.l 
((11/*) = 3Mq), squares for our MGG-9 MC runs ((M*) = 
0.64 Mq). Solid symbols indicate that the model underwent col¬ 
lisional runaway. The models on the left side (low to f values) are 
for MGG-11, the ones on the right for MGG-9. The solid and 
dash-dotted lines indicate where results of GFR04 predict mass- 
segregation-driven core collapse to happen at t = t* = 3 Myr for 
the two different IMFs. 

by GFR04, t C c| rlx < t* = 3 Myr, pre dicts that many more 
of the MGG-11 models considered bv IPortegies Zwart et al] 
should have experienced collisional runaway (see Fig. 1231 . 
This disagreement indicates again that our simple t C c| rlx 
criterion only applies to systems containing a sufficiently 
large number of stars. In particular, for relatively low N*, 
binaries can form through three-body processes and suspend 
core collapse before collisions occur. Furthermore, systems 
of very high concentration and (relatively) low number of 
stars, contain a very low number of high-mass stars in their 
core, initially. Hence, a clean mass-segregation-induced core 
collapse of the type found in our MC simulations may not 
be possible. 

2.5.3 Loss-cone effects 

The MC treatment of collisions is based on the assumption 
that the two neighbouring particles picked up at a given 
step as potential collision partners are members of the same 
(local) distribution. In this approach, one does not actually 
test whether the orbits of the two particles would lead them 
to collide with each other. Instead, it is based on the aver¬ 
age collision probability between such stars in an imaginary 
volume over which the properties of the stellar population 
(density, mass spectrum, velocity distributions...) should 
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Figure 24. Growth of the runaway star. We compare the results 
of simulations using standard collision treatment (with sticky 
sphere approximation) with runs for which the runaway star is 
considered as a fixed central object and the loss -cone effects are 
taken into account following the treatment of iFreitag &; Ben J 
ll2002F) . All runs are for a cluster with Wo = 8, N* = N p = 
1.45 x 10 6 and -Rnb — 0.96 pc. The curves in solid line are for the 
normal collision treatment (K8-14 and K8-15). The dashed curves 
are cases for which a central fixed star of 120 Mq is present from 
the beginning (K8-11, K8-12, K8-17). To obtain the two dash- 
dotted curves, we introduced a central stellar object in a cluster 
which had evolved to deep core collapse using normal collision 
prescription (and no central object). The simulations with normal 
collision treatment were interrupted during the runaway phase 
(models RA095 and RA096). Stellar evolution was not taken into 
account except in one run with loss-cone treatment (K8-17), in 
which case the VMS left the MS at t = 2.8 Myr when its mass 
had reached 3700 Mq. 

be homogeneous. For this method to yield accurate colli¬ 
sion rates, the cluster properties must be relatively constant 
over the range of radii used to estimate the local stellar den¬ 
sity 3 . A more subtle condition, in cases where one particle 
has unique properties and cannot be thought of as a typi¬ 
cal member of some continuous distribution, is that its orbit 
should allow it to explore a significant fraction of this imag¬ 
inary volume (for which the stellar density and collisions 
times are computed) so that its collision probability is also 
relatively homogeneous. 

Unfortunately, these conditions may cease to be fulfilled 
when the runaway star has grown to such a mass that it 
basically stays at the centre of the cluster with little ra¬ 
dial motion in comparison with the size of the orbits of 
nearby stars. The situation then becomes similar to that of 
a massive (or intermediate-mass) black hole at a the cen¬ 
tre of a galactic nucleus (or globular cluster). While the 

3 In the innermost parts of the cluster, this density is computed 
from t he volume of the spherical shell containing 17 nearby par¬ 
ticles jHenonlll973fl . 


VMS grows by merging with stars coming into physical con¬ 
tact with it, the (I)MBH destroys stars venturing into the 
Roche zone and accretes (some of) their material. In both 
cases, stars must come very close to the cluster centre, i.e., 
have very small pericentre distance R pe r i, so the collision or 
disruption rate is dominated by stars on high-eccentricity 
orbits. How these orbits are repopulated by 2-body relax- 
ation is known as the “loss-cone replenishmen t problem” 
(iFrank fc R.eedll97ft iLightman fc Shapirolll977l) . The basic 
difficulty is that relaxation may cause large relative changes 
in .Rperi in only a small fraction of the relaxation time, i.e., 
the time over which the energy of the orbit would be strongly 
affected. To give an illustration, for a Keplerian orbit of ec¬ 
centricity e « 1, The “pericentre relaxation time” is of order 
t rjPer i ~ (1 — e)f r i x where t r i x is the usual (energy) relaxation 
time. Let P or b be the orbital period of the star (from one 
pericentre passage to the next). For regions of phase space 
where P or b <gC f r , P eri, the loss cone is empty, while in the op¬ 
posite regime it is full (since relaxation replenishes loss-cone 
orbits as quickly as stars are removed through interactions 
with the central object). According to standard loss-cone 
theory, in most situations, the disruption rate is dominated 
by stars at the boundary between the two regimes. Hence, 
to include loss-cone physics correctly into stellar dynamical 
models, one should in principle resolve the effects of relax¬ 
ation on timescales much shorter than t r i x , actually as short 
as P or b for stars in the critical regime. The MC scheme used 
here does not allow this directly because the time step St of 
a given particle can only be a function of its present radial 
position R, not of its eccentricity. Furthermore, 5t(R) must 
be an increasing function of R so that setting St = P or b at 
any given R would impose time steps at least as short as this 
for all particles interior to R. To handle the loss-cone prob¬ 
lem in a satisfying way despite of this, an approximate treat¬ 
ment of the effect of relaxation on timescale s shorter than St 
has be en devised; it is described in detail in IFreitag fc Bern] 
( 20021 1. But this approach, as all standard loss-cone theories, 
is based on the assumption that the central object is rigor¬ 
ously fixed at R = 0; in other terms, no account is made of 
its “wandering” through Brownian motion. 

In summary, our present code allows us to follow the 
growth of a VMS using two opposite simplifying assump¬ 
tions. (1) The “standard collisional treatment” which com¬ 
putes collision probabilities with the “nS' co iiV)d” formula 
assuming a small central homogeneous box, ignoring small- 
number effects and the small amplitude of the VMS motion. 
(2) The “loss-cone approach” which neglects this motion al¬ 
together but looks at the possibility of collision between the 
VMS and a given star with (approximate) account of the 
specific orbit of the star and its evolution (due to relax¬ 
ation) on appropriately small timescales. A limitation of the 
second approach, as implemented in the MC code is that it 
only works for interactions with the central object leading 
to disappearance of the star in one pericentre passage. In 
consequence, here, we can only allow mergers but not fly-by 
collisions. 

In Fig. |^] a comparison is made of the results of both 
treatments for a cluster with Wo = 8 , N* = N p = 1.45 x 10 6 
and Rnb = 0.96 pc. Obviously, after the first few collisions, 
the growth of the VMS becomes much slower with the loss- 
cone approach, probably as a result of loss-cone orbits being 
rapidly depleted in the innermost region. In this situation, 
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stellar evolution has a chance to catch up with mergers and 
to terminate the growth of the VMS as the interval of time 
between successive collisions decreases. For this particular 
cluster model, when stellar evolution is allowed, the runaway 
star still grows to more than 3000 Mg before it leaves the 
MS (models K8-16 and K8-17). This happens around t = 
3 Myr but, thanks to collisional rejuvenation, the effective 
age of the VMS (its MS lifetime) is only 1.4 Myr. If the 
collision rate is so reduced after the onset of runaway growth, 
clusters which experience core collapse within just a little 
less than 3 Myr can probably not form a VMS more massive 
than a few hundred Mq before growth is interrupted by the 
evolution of the VMS or normal high-mass stars off the MS. 


The correct VMS growth rate, assuming the same M-R 
relation and collision physics, is very likely between what is 
found with our two different treatments of the central ob¬ 
ject. Which approach is the most accurate is difficult to tell. 
The effects of the wandering of the central object on rates 
of interaction with st ars remain t o be combined with the 
loss-cone theory (see Yound Il977l : iMagorrian fc Tremainel 
Il999l : IsigurdssonlI 200 . 1 - . for short discussions). Considering 
length scales pertaining to the interactions of the VMS 
with the cluster indicates that they take place in a com¬ 
plex regime. We take as an example, simulation K8-29 for 
which we have the best data about the central conditions 
in the runaway phase. The MC simulation indicates that 
the VMS has a wandering radius R w an d of a few 10 -4 
(in A-body units or pc) by the time it has grown heav¬ 
ier than 1000 Mg. This value is in agreement with a sim¬ 
ple analysis based on the assumption of energy equipar- 
tition with the surrounding stars, My MsVwand « (M*)er„, 
with (M*) « 50 Mq, and equal amount of kinetic and po¬ 
tential energy in the cuspy density background of these 
stars, p(R) oc R~ 1 ' 5 . However, even though the motion of 
a massive object in a background of lighter particles in¬ 
teracting with it gravitationally obeys energy equipartition 
approximately l lChatte riee et, al1l2002t iDorband et, a 13 ® 

lLaun fc MerritiKOO^l Mhe energy dissipation introduced by 

collisions can lead to non-equipartition, so the agreement 
may be fortuitous. In any case, the wandering radius of the 
VMS is still much larger than its radius (~ 50 R© ~ 10 -6 pc) 
so it cannot, in principle, be considered strictly fixed at the 
centre. On the other hand, the central number density (es¬ 
timated inside R = 1.2 x 10 ~ 3 for 18 successive snapshots) 
is no ~ 1.5 x 10 s , corresponding to an average distance be¬ 
tween stars of (d) = n^ 1 ^ 3 ~ 2 x 10~ 3 and the influence 
radius of the VMS is formally Rh = GMvms /(3cr 3 ) ~ 1 CU 3 
llFrank fc Reed! Il97fih . These scales are significantly larger 
than J? wan d, pointing to the breakdown of our standard col¬ 
lision treatment already when My ms — 1000 Mq. The two 
approaches we have tried make strong simplifications of a 
situation that can probably only be treated correctly by 
direct orbital integration using, for instance, some hybrid 
numerical code in which the central region of the cluster 
is integrated with a direct A-body scheme, while the outer 
region is modelled in the FP approximation. 


3 SUMMARY AND DISCUSSION 
3.1 Summary of results 

This work is a continuation of our study of the dynami¬ 
cal evolution of young, dense stellar clusters. Our goal is 
to determine under which conditions an IMBH can form 
at the centre of such systems. We are exploring the col¬ 
lisional runaway route. This scenario applies to clusters 
in which the accumulation of massive main-sequence stars 
(M» ~ 10 — 120 Mg) at the centre, through the relaxation- 
driven process of mass segregation, leads to the Spitzer in¬ 
stability, thus triggering core collapse, before these stars 
evolve to become compact objects, i.e., within ~ 3 Myr. As 
the central density of massive stars increases and their ve¬ 
locity dispersion decreases, direct collisions between stars 
must occur, leading to the possibility of successive mergers 
and the runaway formation of one very massive star (VMS, 
M* > 400 Mq). Such a VMS is a possible IMBH progenitor. 

In GFR04, we have established the cluster conditions 
required for the core collapse to be sufficiently rapid. In the 
present work, we have investigated the next stage by imple¬ 
menting stellar collisions. We have followed the core collapse 
and collisional runaway using a Monte Carlo simulation code 
specially devised to account for stellar collisions in a real¬ 
istic way (ME(SSY)**2). In GFR04, 2-body relaxation was 
the only physical process considered in most simulations so 
that, as expected, up to statistical fluctuations, the results 
did not depend on the physical size of the cluster or on the 
number of stars (provided there are enough for the core to 
be resolved). The significant free parameters were the initial 
structure of the cluster (profiles of density, velocity and, in 
some cases, mass segregation) and the IMF. Here, two more 
ingredients, collisions and stellar evolution, are introduced, 
which break this degeneracy. For instance, the ratio of the 
physical central velocity dispersion to the escape velocity 
from a star (hence N,R*/Rh) determines gravitational fo¬ 
cusing and the outcome of collisions and, more importantly, 
the ratio tcoii/trix- 

We have computed more than 100 models, varying the 
cluster size (7?h in the range 0.03 — 5pc), mass (A* = 
10 5 — 10 8 ) and concentration. For the latter parameter, we 
have considered King models with Wo = 3 and Wo = 8 . 
For most simulations we have used a Salpeter IMF extend¬ 
ing from 0.2 to 120 Mq. We treated collisions either in the 
sticky-sphere approximation or by making use of prescrip¬ 
tions for merger conditions and mass loss derived from a 
large set of SPH simulations. We used up to 9 x 10 6 par¬ 
ticles (one simulation) but, because of the large number of 
simulations and the relatively large amount of CPU time 
required to track the runaway growth (requiring very short 
time steps), we had to perform most high-A* runs using a 
number of particles lower than A*. Taking advantage of the 
statistical nature of the MC algorithm, our code scales phys¬ 
ical processes so that each particle may represent a given 
(fixed) number of stars, not necessarily one. This feature 
may legitimately be questioned when a runaway happens as 
one expects only one star to experience it. However, by vary¬ 
ing the number of particles (but keeping A* fixed) for a few 
cluster simulations, we have checked that using A p < A* 
does not appear to produce spurious results (such as, e.g., 
preventing the runaway). 

We found that runway VMS formation occurs in all 
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cases for which the core collapse time is shorter than the 
MS lifetime of massive stars, as predicted in GFR04. Fur¬ 
thermore, for very massive clusters (> 10' Mg), the core 
collapse is facilitated (and sometimes driven) by collisions 
themselves, an effect which extends the runaway domain to 
clusters of larger size than predicted by an analysis based 
solely on relaxation. Such systems have velocity dispersions 
in excess of 300 km s^ 1 but we find that the collisional mass 
loss and the reduced merger cross section cannot prevent 
runaway, even at ~ 1000 km s , 

We never observe the formation of more than one VMS. 
Only in clusters that are initially collisional do collisions pro¬ 
duce a high-mass tail in the stellar mass spectrum. Most of 
the mass eventually incorporated into the VMS originates 
in stars near the top of the IMF. Typically, for clusters that 
are not initially collisional, > 90 % of the VMS mass is con¬ 
tributed by stars in the 60 — 120 Mg range. This is due to the 
larger cross section of massive stars and, more importantly, 
to the fact that collisions take place in the collapsed core, 
which is dominated by these stars. In such clusters, most 
stars that experience a collision do so only once, namely 
when they merge with the growing VMS. In very dense clus¬ 
ters where collisions are happening from the beginning of the 
evolution, the contribution of lighter stars is slightly higher 
and a significant fraction of stars that merge with the VMS 
have experienced at least one previous collision. 

Changing the parameters of the IMF modifies the con¬ 
dition for rapid core collapse but has little qualitative effect 
on the runaway pro cess itself . In partic ular, clusters with 
a field Kroupa IMF llKroupa et al.lrl993ll . which has an ex¬ 
ponent a = 2.7 for massive stars, do experience runaway 
in essentially the same way as those with a Salpeter IMF 
(a = 2.35), despite a much smaller total number of massive 
stars. If the upper cutoff of the IMF is significantly lower 
than the assumed M max = 120 Mg, core collapse will take 
longer but for M max < 30 Mg, this is over-compensated by 
a longer MS lifetime. From one simulation we performed for 
M max = 10 Mg, we predict that the region of parameter 
space leading to runaway is similar to what we found for 
Af max = 120 Mg. 

We found that in most circumstances, once the run¬ 
away has started, the average time between successive merg¬ 
ers becomes much shorter than the MS thermal relaxation 
timescale for the VMS, Ikh- Hence, its structure and re¬ 
sponse to further collisions may be greatly affected by its 
prior collision history. Modelling the complex interplay be¬ 
tween stellar dynamics, collisions, stellar structure and stel¬ 
lar evolution which should eventually determine the final 
mass of the VMS is presently not possible. Among the ef¬ 
fects with potentially strong bearing on the runaway pro¬ 
cess is the depletion of loss-cone orbits that may come into 
play before the VMS has reached ~ 1000 Mg as, being more 
massive than surrounding stars, it becomes confined to a 
very small central region of the cluster. Assuming a strictly 
fixed central object and using the approximate treatment 
of loss-cone physics built into our MC code (for the study 
of massive black holes in galactic nuclei), we have found a 
significant reduction in the VMS growth rate (after the first 
2 — 4 mergers) but, in the few cases considered, the VMS still 
grew to a few thousand Mg before evolving off the MS and 
the average time between collisions, although 3 — 10 times 
longer, remained below Ikh most of the time. 


Our simulations allow us to establish the conditions for 
the onset of runaway collisions and to study the early VMS 
growth. If globular clusters are formed with a concentra¬ 
tion equal or higher than a Wo = 8 King model (and if 
primordial binaries cannot prevent the runaway process), a 
significant fraction (> 20 %, see Fig. 1 of Paper I) may ex¬ 
perience this phase in the first few million years of their 
dynamical evolution. Typically, a cluster containing ~ 10 6 
stars at this concentration needs to have a half-mass radius 
smaller than i?h ~ 2 pc to experience runaway, but a proto- 
galactic nucleus with 10 s stars has to be more compact than 
R h ~ 0.8 pc. 

3.2 Discussion 

Many aspects of the runaway collision scenario and how we 
have approached it in our numerical simulations have been 
discussed already. In particular, the choice of physical ingre¬ 
dients considered and the uncertainties and simplifications 
involved in their treatment have been presented in Paper I. 
In the results section of the present paper, we have discussed 
other important uncertain aspects of the process, namely the 
impact of collisions on the structure and growth of the VMS 
and the possible effects of loss-cone depletion. In this last 
subsection, we consider the role of interstellar gas. 

3.2.1 Role of residual gas in young clusters 

Because the runaway scenario presented here has to oper¬ 
ate within the first few Myr of the cluster life, it has, in 
fact, to be considered in the framework of cluster formation, 
which takes place on a similar timescale. No only does it 
mean, as mentioned above, that stars below 2 — 3 Mg may 
still be on the pre-MS, with much larger radii than adopted 
here, but, more importantly, that a significant amount of 
residual gas may be present. Observations show that when 
a cluster forms, no t more than 30 % of the gas is eventually 
turned into stars llLadalll999lf but that clusters like R136 
or NGC 3603, as young as 1 — 2 Myr, are already devoid of 
gas fe.g.. IStolte et alJl2004 '). In such clusters, the expulsion 
of gas (originally the dominant mass component) is driven 
by the ionising radiation and winds of OB stars and occurs 
impulsively. A large fraction of the stars are then lost from 
the cluster and the remaining bound cluster swells quit e dra¬ 
matically llKroupa et all200lHBoilv' fc Kroupai2003al O) , an 
event likely to terminate the collisional phase after ~ 1 Myr 
rather than ~ 3 Myr. In clusters with an escape veloc¬ 
ity much higher than the typical expansion velocity of the 
ionised gas, i.e., ~ lOkms , complete expulsion of the gas 
probably only occurs when the first SN explodes. 

The faster evolution due to the much smaller size of 
the cluster in the embedded phase counterbalances the pos¬ 
sible reduction in the time available for collisional run¬ 
away. Assuming that, when gas is expelled, the mass of 
the cluster M c \ (including gas) is reduced by a factor of 
~ 4 and its size Rh expands ~ 10 times, the true initial 
relaxation time (t r i x oc crj)n _1 oc M^{ 2 R^ 2 ) should be of 
order 4 times shorter than what is inferred from observa¬ 
tions of the gas-free cluster. Furthermore, the collision time 
(tcoii oc <7„n _1 oc M~^ 2 Ry/ 2 ) may have been some 1000 
times shorter, suggesting that a large fraction of clusters 
may actually be collisional at birth. 
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The standard theory of star formation, which assumes 
spherically symmetric accretion of gas on t o the proto¬ 
star, fails for stars more massive than 10 Mq (IYorkell2004 
iBallv fc Zinneckerll2005f ). This problem led to the proposal 
that the truly primordial IMF may not extend past 10 Mq 
and that mergers may be —at least pa rtially— responsible 
for the formation of more massive stars l| Bonnell et al ] 1991 
iBonnell fc Batell2002l : IBallv fc Zinneckerll2005f) . This would 
also naturally explain why these massive stars are preferen¬ 
tially found in the central regions of clusters that are pos¬ 
sibly too young to have experienced relaxational mass seg¬ 
regation feonnell fc Daviea 11994 Ide Griis et al J l2002fl . To 
occur on a sufficiently short timescale, collisions require cen¬ 
tral stellar densities orders of magnitude higher than what 
is observed in those young clusters but similar to what it 
may have been in the embedded phase. 

Realistic simulations of embedded clusters should in¬ 
clude the role of gas, since its mass is likely to be more 
than twice that in stars. Especially important are the in¬ 
creased velocity dispersion and accretion from the ambient 
gas. The pioneering SPH/N-body simulations of very early 
cluste r evolution could only handle up to a fe w hundred 
stars feonnell et alJ 120041 : feate fc Bonnelll 1200511 . Our MC 
results (run K3-20 presented in Section [2.5.111 suggest that, 
in much larger systems, collisions would lead to the forma¬ 
tion of just one (very) massive star rather than a continuous 
spectrum spreading from 10 Mq to ~ 100 Mq . However, this 
result applies to clusters driven to core collapse by mass seg¬ 
regation. When the cluster is initially collisional, more stars 
should take part in collisions and a high-mass spectrum may 
develop in a more orderly way (see Fig. Il.it . Also, the colli¬ 
sional phase may be terminated before a runaway collision 
sequence starts through gas expulsion by massive stars. 


3.2.2 Gas retention in galactic nuclei 

Our predictable result that evolution of the massive stars 
off the MS terminates the evolution towards core collapse 
stems from our assumption that all gas released by the stars 
is completely (and instantaneously) lost from the cluster. 
If the gas is not totally expelled, there may be no core re¬ 
expansion. The core may instead remain very dense. This 
opens the possibility for a longer-term collisional evolution 
but now there will be plenty of stellar BHs around that may 
destroy MS stars. 

Over a long timescale, it is possible t hat star forma¬ 
tion w ill convert the ga s into new stars t Sanders! 197(1 : 
IQuinlan fc Shapirolll990 l. All this becomes obviously quite 
intricate but would have to be included to treat realisti¬ 
cally proto-galactic nuclei. Nowadays, multi-phase dynam¬ 
ical simulations, including interstellar gas (possibly in two 
phases), stars and dark matter, and the matter and energy 
exchanges between these components, are commonly per¬ 
formed to investigate the formation and evolution of gala x¬ 
ies (e.g., iMihosJ^H emguisj 1199(1 lAbadi et, alJ l2003al lbt 
ISpringel fc Hernnuis tEM . These works are based on the 
use of collisionless IV-body codes coupled with SPH. In¬ 
cluding detailed gas dynamics in models of relaxing clus¬ 
ters is extremely chalenging because the gas should pos¬ 
sess a complex, non spherically symmetri c structure, evolv- 
ing on timescales much shorter than t r i x llCoker et al.lil9~9 9: 


[Williams et al.lll99fll : iRockefeller et, all 120041 : ICuadra et alJ 

Note that, in clusters with very high velocity disper¬ 
sions, where collisions may significantly contribute to gas 
production, some of it will probably be retained. Indeed 
this can be seen in the extreme situation of a completely 
disru ptive collision (a rare occurrence; see [Freitag fc Bend 
12005th In this case, the total energy of the colliding stars 
is their kinetic energy at (relatively) large separation minus 
the sum of their individual binding energy (a positive quan¬ 
tity, by definition). None of the stars can have a velocity 
higher than the cluster’s escape velocity 14 sc and because 
the binding energy is positive, the energy available per unit 
mass of stellar gas has to be smaller than 1/2 V BSC so at least 
some gas must stay in the cluster to conserve energy. 
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Table 1. Properties of Simulated Clusters 


Name 

W 0 

IV, 

N p 

-Rnb 

(pc) 

Coll. 

Peculiarities 

tra 

C t fp ) 

tr& 

(Myr) 

K3-1 

3 

to 5 

10 s 

0.3 

ss 

1 - 120 Mg 

4.38 

X 

10 -2 

2.77 

K3-2 

3 

to 5 

10 s 

0.4 

ss 


1.08 

X 

10 -2 

2.25 

K3-3 

3 

10 s 

10 s 

0.4 

ss 

1 - 120 Mg 





K3-4 

3 

10 s 

10 s 

0.5 

ss 






K3-5 

3 

to 5 

10 s 

0.5 

ss 


1.05 

X 

10“ 2 

3.03 

K3-6 

3 

10 s 

10 s 

0.5 

ss 


1.07 

X 

10 -2 

3.08 

K3-7 

3 

10 s 

10 s 

0.5 

ss 


1.02 

X 

10" 2 

2.94 

K3-8 

3 

to 5 

10 s 

0.5 

ss 

1 - 120 Mg 





K3-9 

3 

3 X 10 5 

3 x 10 5 

0.03 

ss 


9.89 

X 

10 -3 

0.0634 

K3-10 

3 

3 x 10 5 

3 x 10 5 

0.1 

ss 


1.06 

X 

10 -2 

0.414 

K3-11 

3 

3 x 10 5 

3 x 10 5 

0.1 

SPH 


1.04 

X 

10~ 2 

0.406 

K3-12 

3 

3 x 10 5 

3 x 10 5 

0.2 

SS 


1.05 

X 

10" 2 

1.16 

K3-13 

3 

3 X 10 5 

3 x 10 5 

0.3 

ss 


9.60 

X 

10 -3 

1.95 

K3-14 

3 

3 X 10 5 

3 x 10 5 

0.3 

ss 

Constant 5 % mass loss 

1.03 

X 

10 -2 

2.09 

K3-15 

3 

3 x 10 5 

3 x 10 5 

0.3 

ss 

Constant 10 % mass loss 

1.03 

X 

10" 2 

2.09 

K3-16 

3 

3 x 10 5 

3 x 10 5 

0.3 

ss 

Constant 30 % mass loss 





K3-17 

3 

3 X 10 5 

3 x 10 5 

0.3 

ss 

Fixed central object 

1.01 

X 

10 -2 

2.05 

K3-18 

3 

3 x 10 5 

3 x 10 5 

0.3 

X 

5 % mass loss 

9.35 

X 

10 -3 

1.9 

K3-19 

3 

3 x 10 5 

3 x 10 5 

0.3 

ss 

Constant R * for VMS 

9.66 

X 

10" 3 

1.96 

K3-20 

3 

3 x 10 5 

3 x 10 5 

0.3 

ss 

0.2 - 10 Mg 

9.04 

X 

10 -2 

19.9 

K3-21 

3 

3 X 10 5 

3 x 10 5 

0.3 

ss 

Kroupa 0.08 — 120 Mq 

1.02 

X 

10 -2 

2.54 

K3-22 

3 

3 X 10 5 

3 x 10 5 

0.3 

SPH 


8.88 

X 

10 -3 

1.8 

K3-23 

3 

3 x 10 5 

3 x 10 5 

0.4 

SS 


9.50 

X 

10 -3 

2.97 

K3-24 

3 

3 x 10 5 

3 x 10 5 

0.4 

SPH 


9.28 

X 

10" 3 

2.9 

K3-25 

3 

3 X 10 5 

3 x 10 5 

0.4 

ss 


9.69 

X 

10~ 3 

3.02 

K3-26 

3 

3 X 10 5 

3 x 10 5 

0.4 

ss 


1.02 

X 

10 -2 

3.18 

K3-27 

3 

3 x 10 5 

3 x 10 5 

0.4 

ss 

Check for orbit overlap 

9.22 

X 

10 -3 

2.85 

K3-28 

3 

3 x 10 5 

3 x 10 5 

0.5 

ss 






K3-29 

3 

3 X 10 5 

3 x 10 5 

0.5 

ss 

0.2 - 10 Mg 





K3-30 

3 

10 6 

10 6 

0.1 

ss 

Kroupa 0.01 — 120 M© 

3.50 

X 

10 -3 

0.327 

K3-31 

3 

10 6 

10 6 

0.2 

ss 


7.07 

X 

10" 3 

1.87 

K3-32 

3 

10 6 

10 6 

0.2 

SPH 


8.00 

X 

10" 3 

1.39 

K3-33 

3 

10 6 

10 6 

0.3 

ss 


7.94 

X 

10 -3 

2.54 

K3-34 

3 

10 6 

10 6 

0.4 

ss 






K3-35 

3 

3 x 10 6 

3 x 10 5 

0.3 

ss 






K3-36 

3 

3 x 10 6 

3 x 10 6 

0.1 

ss 


6.44 

X 

10" 3 

0.617 

K3-37 

3 

3 X 10 6 

3 x 10 6 

0.2 

ss 


6.84 

X 

10 -3 

1.85 

K3-38 

3 

3 x 10 6 

3 x 10 6 

0.2 

SPH 


7.25 

X 

10 -3 

1.97 

K3-39 

3 

3 x 10 6 

3 x 10 6 

0.3 

SS 






K3-40 

3 

10 7 

3 x 10 5 

0.1 

ss 


7.75 

X 

10 -3 

1.21 

K3-41 

3 

10 7 

3 x 10 5 

0.1 

ss 


7.53 

X 

10 -3 

1.18 

K3-42 

3 

10 7 

3 x 10 5 

0.2 

ss 


8.79 

X 

10 -3 

3.89 

K3-43 

3 

10 7 

3 x 10 5 

0.25 

ss 






K3-44 

3 

to 7 

3 x 10 5 

0.3 

ss 






K3-45 

3 

10 7 

3 x 10 5 

0.4 

ss 






K3-46 

3 

10 7 

10 6 

0.1 

SPH 


6.41 

X 

10 -3 

0.998 

K3-47 

3 

10 7 

10 6 

0.2 

ss 


6.27 

X 

10 -3 

2.76 

K3-48 

3 

10 7 

10 6 

0.25 

ss 






K3-49 

3 

3 x 10 7 

3 x 10 5 

0.1 

ss 


5.24 

X 

10" 3 

1.3 

K3-50 

3 

10 8 

3 x 10 5 

0.03 

ss 


3.49 

X 

10 -4 

0.0237 

K3-51 

3 

10 8 

3 x 10 5 

0.04 

SPH 


2.99 

X 

10 -3 

0.312 

K3-52 

3 

10 8 

3 x 10 5 

0.05 

SS 


7.36 

X 

10~ 4 

0.107 

K3-53 

3 

10 8 

3 x 10 5 

0.05 

SS 

No relaxation 

8.90 

X 

10 -4 

0.13 

K3-54 

3 

10 8 

3 x 10 5 

0.1 

ss 


1.94 

X 

10 -3 

0.801 

K3-55 

3 

10 8 

3 x 10 5 

0.2 

ss 


3.41 

X 

10 -3 

3.98 

K3-56 

3 

10 8 

3 x 10 5 

0.2 

SPH 






K3-57 

3 

10 8 

3 x 10 5 

0.25 

SS 
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Table 1—Continued 


Name 

W 0 

iV* 

N p 

Rnb 

(pc) 

Coll. 

Peculiarities 

tra 

( T FP) 

fra 

(Myr) 

K3-58 

3 

10 s 

3 x 10 s 

0.3 

SS 






K3-59 

3 

10 8 

3 x 10 5 

0.4 

ss 






K3-60 

3 

10 s 

3 x 10 5 

0.5 

SS 






K3-61 

3 

10 s 

10® 

0.04 

ss 


5.05 

X 

10 -4 

0.0524 

K3-62 

3 

10 8 

10® 

0.1 

SPH 


3.19 

X 

10" 3 

1.31 

K3-63 

3 

10 8 

10® 

0.1 

SPH 


3.20 

X 

10" 3 

1.31 

K3-64 

3 

10 8 

10® 

0.1 

SS 


1.90 

X 

10 -3 

0.78 

K3-65 

3 

10 8 

10® 

0.2 

SPH 






K8-1 

8 

3 x 10 s 

3 x 10 5 

0.9 

SS 


3.97 

X 

10 -4 

0.418 

K8-2 

8 

3 x 10 s 

3 x 10 5 

2 

SS 


4.63 

X 

10 -4 

1.62 

K8-3 

8 

3 x 10 s 

3 x 10 5 

3 

SS 






K8-4 

8 

3 x 10 5 

3 x 10 5 

4 

ss 






K8-5 

8 

3 x 10 5 

3 x 10® 

1.2 

SPH 


3.89 

X 

10- 4 

0.618 

K8-6 

8 

3 x 10 s 

3 x 10® 

4.7 

SPH 






K8-7 

8 

10® 

10® 

0.1 

SS 


4.00 

X 

lO" 4 

0.0246 

K8-8 

8 

10® 

10® 

1 

SS 


4.32 

X 

10 -4 

0.841 

K8-9 

8 

10® 

10® 

2 

SS 


4.82 

X 

lO" 4 

2.65 

K8-10 

8 

10® 

10® 

3 

SS 






K8-11 

8 

1.4 x 10® 

1.4 x 10® 

0.96 

SS 

Fixed central object 

4.02 

X 

10 -4 

0.843 

K8-12 

8 

1.4 x 10® 

1.4 x 10® 

0.96 

SS 

Fixed central object 

3.91 

X 

10 -4 

0.819 

K8-13 

8 

1.4 x 10® 

1.4 x 10® 

0.96 

SS 

Fixed central object, tid. trunc. 

3.75 

X 

lO" 4 

0.786 

K8-14 

8 

1.4 x 10® 

1.4 x 10® 

0.96 

SS 


4.66 

X 

10 -4 

0.976 

K8-15 

8 

1.4 x 10® 

1.4 x 10® 

0.96 

SS 


3.98 

X 

io- 4 

0.834 

K8-16 

8 

1.4 x 10® 

1.4 x 10® 

0.96 

SS 

Fixed central object, stell. evol. 

4.19 

X 

10 -4 

0.878 

K8-17 

8 

1.4 x 10® 

1.4 x 10® 

0.96 

SS 

Fixed central object, stell. evol. 

3.91 

X 

10- 4 

0.819 

K8-18 

8 

1.4 x 10® 

1.4 x 10® 

0.96 

SS 


4.32 

X 

10 -4 

0.904 

K8-19 

8 

3 x 10® 

3 x 10® 

0.2 

SPH 


3.31 

X 

10- 4 

0.0892 

K8-20 

8 

3 x 10® 

3 x 10® 

1 

SS 


3.43 

X 

10 -4 

1.03 

K8-21 

8 

3 x 10® 

3 x 10® 

1 

SPH 


3.43 

X 

10- 4 

1.03 

K8-22 

8 

3 x 10® 

3 x 10® 

1 

SS 


3.70 

X 

10" 4 

1.12 

K8-23 

8 

3 x 10® 

3 x 10® 

1 

SS 


3.74 

X 

10 -4 

1.13 

K8-24 

8 

3 x 10® 

3 x 10® 

1 

SS 

Fixed central object 

2.99 

X 

10 -4 

0.9 

K8-25 

8 

3 x 10® 

3 x 10® 

1.2 

SPH 


4.22 

X 

10 -4 

1.65 

K8-26 

8 

3 x 10® 

3 x 10® 

1.4 

SS 


4.27 

X 

10 -4 

2.13 

K8-27 

8 

3 x 10® 

3 x 10® 

2 

SS 






K8-28 

8 

3 x 10® 

3 x 10® 

4.7 

SPH 






K8-29 

8 

9 x 10® 

9 x 10® 

1 

SS 


3.78 

X 

10 -4 

1.78 

K8-30 

8 

10 7 

3 x 10® 

1 

SS 


3.61 

X 

10 -4 

1.78 

K8-31 

8 

10 7 

3 x 10® 

1.4 

SS 


3.12 

X 

10 -4 

2.55 

K8-32 

8 

10 7 

3 x 10® 

1.4 

SS 


3.86 

X 

10 -4 

3.15 

K8-33 

8 

10 7 

3 x 10® 

2 

SS 






K8-34 

8 

3 x 10 7 

3 x 10® 

1 

SS 


3.14 

X 

10 -4 

2.45 

K8-35 

8 

3 x 10 7 

3 x 10® 

1.4 

SS 






K8-36 

8 

10 8 

3 x 10® 

0.1 

SS 


7.49 

X 

10 -5 

0.0308 

K8-37 

8 

10 8 

3 x 10® 

0.6 

SS 


2.40 

X 

10 -4 

1.45 

K8-38 

8 

10 8 

3 x 10® 

0.8 

SS 


2.74 

X 

10' 4 

2.55 

K8-39 

8 

10 8 

3 x 10® 

1 

SS 






K8-40 

8 

10 8 

3 x 10® 

1.4 

SS 






K8-41 

8 

10 8 

3 x 10® 

2 

SS 






MGG9-K8 

8 

2.7 x 10® 

2.7 x 10® 

2.6 

SPH 

Model for MGG-9 





MGG9-K9 

9 

2.7 x 10® 

2.7 x 10® 

2.6 

SPH 

Model for MGG-9 

1.55 

X 

10' 4 

1.95 

MGG9-K12a 

12 

2.7 x 10® 

2.7 x 10® 

2.6 

SPH 

Model for MGG-9 

5.10 

X 

10-® 

0.0644 

MGG9-K12b 

12 

2.7 x 10® 

2.7 x 10® 

2.6 

SPH 

Model for MGG-9 

5.98 

X 

10-® 

0.0755 


Note. — When not otherwise mentioned in the ”Peculiarities” column, clusters were modeled with a Salpeter (a = 2.35) 
IMF extending from 0.2 to 120 Mq, our standard M-R relation and prescription for collisional rejuvenation and MS lifetime. 
The data listed in the column entitled ”Coll.” indicates the type of collision treatment used. ”SS” stands for sticky-sphere 
approximation; ”SPH” for the prescriptions for mass loss and merger derived from the SPH simulations of Freitag &; Benz 
(2005); ”X” is a special case for which a simple merger criterion was used (A me rg = —1.14 — 0.5(Z V —3.0), see Section 3.2.5 of 
Paper I) and a flat mass-loss rate of 5 % was assumed in all cases. t ra is the time at which runaway growth of a VMS started. 
It is given in Fokker-Planck time units (Tpp) and Myr. No value is given for clusters which didn’t experience collisional 
runaway. 











